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The fractional symmetric rigid rotor 
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Abstract. Based on the Riemann fractional derivative the Casimir operators and multiplets for the frac- 
tional extension of the rotation group SO(n) are calculated algebraically. The spectrum of the correspond- 
ing fractional symmetric rigid rotor is discussed. It is shown, that the rotational, vibrational and 7-unstable 
limits of the standard geometric collective models are particular limits of this spectrum. A comparison with 
the ground state band spectra of nuclei shows an agreement with experimental data better than 2%. The 
derived results indicate, that the fractional symmetric rigid rotor is an appropriate tool for a description 
of low energy nuclear excitations. 

PACS. 21.60.Fw Group theory nuclear physics - 21.10.k Nuclear energy levels 



1 Introduction 

The concept of fractional calculus has inspired mathemati- 

■ cians since the days of Leibniz [?]- [2]. 

In physics, early attempts to use this concept in the 
field of applications had been studies on non-local dy- 
namics, e.g. anomalous diffusion or fractional Brownian 
motion [3], [4]. 

During the last decade, remarkable progress has been 
made in the theory of fractional wave equations[5]-[5]. But 
until now, there exists no application of the fractional cal- 

■ cuius in nuclear theory. 

We therefore will first introduce the necessary tools 

■ for a fractional extension of the standard definition of an- 
gular momentum operators, derive the Casimir operators 
and classify the multiplets of a fractional extension of the 
standard rotation group SO(n). 

We will then discuss the properties of the correspond- 
ing fractional symmetric rigid rotor spectrum and the re- 
sults of its application to the low energy excitation ground 
state band spectra of even-even nuclei. 



2 Fractional quantum mechanical observables 

The transition from classical mechanics to quantum me- 
chanics may be performed by canonical quantisation [S], 
[TO] . The classical canonically conjugated observables x 
and p are replaced by quantum mechanical observables X 
and P, which are introduced as derivative operators on a 
Hilbert space of square integrable wave functions /. The 
space coordinate representations of these operators are 



Xf(x) = xf(x) 



(1) 



Pf(x) = -ihd x f(x) 



where X and P fulfill the commutation relation 



X,P 



= ih 



(2) 



(3) 



We will now describe a generalisation of these operators 
from integer order derivative to arbitrary order derivative. 

The definition of the fractional order derivative is not 
unique, several definitions e.g. the Caputo, Weyl, Riesz, 
Griinwald |llj-[14j fractional derivative definition coex- 
ist. For our purpose, we apply as a specific representa- 
tion of the fractional derivative operator the widely used 
Riemann-Liouville fractional derivative [15j of order a, 
where a is a real number in the interval [0, 1]: 



D%(x)f(x) 



1 



d 



r(l - a) dx j 



dt(x-S)- a f@ (4) 



Since ((4]) is valid for x > only, we define the variable \ 
and the corresponding fractional derivative as: 



sign(x) \x\ a 
sign(x) D%(\x\) 



(5) 
(G) 



With the definitions (U])-© we are able to calculate frac- 
tional derivatives for real numbers R, e.g. for f(x) — x" 
we obtain: 



rn 



X 



r(l + (1/ - l)a) 



X 



va > — 1 



(7) 



Furthermore, we are able to define Riemann- Taylor series 
fix) on R: 
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fix) = \x\ 



n=0 



(8) 
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To construct a Hilbert space with these functions /(x) we 
use the Riemann-Liouvillc fractional integral 15J of order 
a: 

4 a (*)/(*) = if^(i- a ~V(0 (9) 

Just like the Riemann-Liouville fractional derivative, this 
integral is valid for x > only. Therefore, we define an 
integral operator on R: 



du = sign(x)I"(\x\) 



The fractional scalar product is then given as: 



< /l<7 >= / 



duf*(x)g(x) 



(10) 



(11) 



With these definitions, the fractional quantum mechani- 
cal observables X and P finally read in space coordinate 
representation as: 



xf(x) 



Pf(x) 



(l-a) 



Xf(x) (12) 
mcDzfix) (13) 



The attached factors (li/mc)' 1 " ' and (h/mc) a mc respec- 
tively ensure correct length and momentum units. For 
a = 1 these definitions reduce to (pj and ([2]). 

Since the Leibniz' product rule is not valid in its origi- 
nal form anymore for fractional derivatives [16) , commuta- 
tion relations depend on the specific choice of the function 
set, they are acting on. For a function set {/(x) = X v }> 
using ([7]), the fractional conjugated operators satisfy the 
following commutation relation: 



X,P 



= —ih 



r(i 



r(l + (i/-l)a) 
ihc(v, a) 



r(l + (u + l)a) i 
r{l + va) 



(J4) 
(15) 



In (|15p we have introduced —ihc(v,a) as a short hand 
notation of the commutator. For a — 1 this commutator 
is equivalent with ©. 

The results derived so far may easily be extended to 
the multi dimensional case. For the Euclidean space of N 
particles we obtain the following sets of fractional conju- 
gated operators in space coordinate representation: 



{Xi} = { i — i 



mc J 



Xi} 



{Pi} = {-i — mcDi} 
mc , 



(16) 

t = l,...,30/7) 



where we have introduced the abbreviation Di for . 

On a function set of the form {/(xij X2, X3Jv) = 
rii^Xr}' ^ ne fractional conjugated operators satisfy the 



following commutation relations: 

Xi , Xj = 
Pi, Pi] =0 

.V,./' ,1,6 J r (l + vi<*) 
- ' J \ 3 \r{\ + { Vi -l)a) 

= -ihSi 3 c(vi, a) 



(18) 
(19) 



r(i + (uj + i 

r(l + via) 



(21) 



With definitions (jT6J) , (|T7J) , we are now able to quantise 
classical operators, e.g. the classical angular momentum 
^classical _ _ uge Q r canonical quantisa- 

tion procedure [5], [ID]. 

3 Classification of angular momentum 
eigenstates 

According to the results from the previous section, we 
define the generators of infinitesimal rotations in the 
planc (i,j = 1, ...,3-ZV), with N being the number of par- 
ticles): 



= -ifi (xiDj - Xi A) 

On a function set 

3N 

{/(xi,x2, -,x3jv) =n^ 



(22) 



(23) 



the commutation relations for Lij are isomorph to an ex- 
tended fractional SO a (3N) algebra: 

\Lij^Lrnn\ — iTlc{y,Ct) {^&iml j jn~\~&jnl j irn Sin-Ljm $jmLin) 

(24) 

Consequently, we can proceed in a standard way [18j,[19]. 
We define a set of Casimir operators A\ , where the index 
k indicates the Casimir operator associated with SO a (k): 



4 = jE(^) 2 



fc = 2,...,3iV (25) 



1,3 



which indeed fulfills the relations [A 2 N , L^] = and suc- 
cessively [A 2 k , A 2 k ,] = 0. Therefore the multiplets of SO a (3N) 
are classified according to the group chain 

SO a {3N) D SO a (3N - 1) D . . . D SO a (3) D SO a (2) 

We use Einstein's summmation convention and introduce 
the metric of the Euclidean space to be gij = 5ij for rais- 
ing and lowering indices. The explicit form of the Casimir 
operators is then given by 

A\ = +X i XiP j P j - ihc{v, a) (k - l)P Pj - X'PP.Pj 

i = l,...,k (27) 
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The classical homogeneous Euler operator is defined as 
xidxi + x 2dx 2 + ---~\~ x kdx k ■ We introduce a generalisation of 
the classical homogeneous Euler operator J k for fractional 
derivative operators 



satisfy the Laplace equation D l Dif = and are nor- 
malised in the interval [—1, 1]. 

We propose the quantisation condition: 

v = n n = 0, 1, 2, ... 



(36) 



1. 



(28) 



Where n is a non-negative integer. This specific choice 

reduces to the classical quantisation condition for the case 
With the generalised homogeneous Euler operator the Casimir^ _ ^ 

operators are: Q n tnig jjilbert space H a , the generalised homoge- 

neous Euler operator J k is diagonal and has the eigen- 



A\ = ■ XX.I'T, 



(c{v,a) (k~2)r k +r k r k 



i — 1, k values lk(a,n) 
(29) 



From this equation it follows, that the Casimir operator 
is diagonal on a function set /, if the generalised homoge- 
neous Euler operator is diagonal on / and if / fulfills the 
Laplace-equation D l Dif — 0. 

We will show, that the generalised homogeneous Euler 
operator is diagonal, if / fulfills an extended fractional ho- 
mogeneity condition, which we will derive in the following. 
For that purpose, we first verify, that the generalised ho- 
mogeneous Euler operator is diagonal on a homogeneous 
function in the one dimensional case. 

We apply the following scaling property, which is valid 
for the fractional derivative [IB], with A = sign(X) \X\ a : 

DxfCm = XD- X ^(\x) (30) 
Homogeneity of a function in one dimension implies: 



/(Ax) = A"/(X) 



VOL > — 1 



(31) 



Applying the derivative operator on (f3"Tj) . using ([3T)|) 
and ((Tj) leads to: 



r(l + va) 



r(l + {v-l)a) 
Setting X = 1 reduces (|32|) to: 

r(i + vol) 



xD^f(x) 



r{l + (v-l)a) 



fix) 



(33) 



The left hand side of (|3"3"|) is nothing else but the gener- 
alised homogeneous Euler operator (|28]) in one dimension 
J e 

Therefore, for the multi dimensional case, on a function 
set of the form 

{/(Xi,X2,--,Xfc) =^auiv 2 ...u k XiX2 2 -Xk h } (34) 



the Euler operator J k is diagonal, if the Vi fulfill the fol- 
lowing condition: 



E 



r(i 



r(i 



va) 



-r(l + (v t -l)a) r(l + (v-l)a) 



(35) 



This is the fractional homogeneity condition. 

Hence we define the Hilbert space H a of all functions 
/, which fulfill the fractional homogeneity condition (|35p . 



l k (a,n) 



n = 0,1,2, 



(37) 



r(l + (n-l)a) 
The eigenvalues of the Casimir-operators on H a follow as: 
A 2 f = hl 2 (a,n)f 

)(z fc (o,n) + c(n,a)(fc-2))/ fc = 3,. 



A k f = h lk(a,n 
with 

h(ot,n) > lk-i(a,n) > 



A 2 



^12 



j\ z it follows 



J z (a) \JM> 



r(l + Ma) 



\JM> 



J 2 (a) | JM> = h 



r(l + (J-l)a) 
J = 0,+l. 



\JM> 
+2,... 



(38) 
,(39^ 



>| ±l 2 (a,n) |>0 (40) 



We have derived an analytic expression for the full spec- 
trum of the Casimir operators for the SO a (3N), which 
may be interpreted as a projection of function set J34|) on 
to function set ([23]) determined by the fractional homo- 
geneity condition (f35l) . 

For the special case of only one particle (TV = 1), we 
can introduce the quantum numbers J and M, which now 
denote the J-th and M-th eigenvalue h(a, J) and h{oL, M) 
of the generalised homogeneous Euler operators Jf and 
Jf respectively. The eigenfunctions are fully determined 
by these two quantum numbers / =| JM >. 

With the definitions J z (a) = A% = L\% and J 2 (a) 

^13 + ^2 



(M - l)a) 
M = 0,±1,±2,...,±J (41) 
2 r(l + (J + l)a) 



(42) 



For a — 1 equations (|4Tj) and (|42|) reduce to 

1) \JM> = HM\JM> M = 0,±l,±2,...(i3| 



J z {a 
J 2 (a 



1) \JM> = h 2 J(J + 1) \JM> J = 0,+l,+^44) 



and J 2 (a = 1), J z (a = 1) reduce to the definitions of 
J 2 ,J Z used in standard quantum mechanical angular mo- 
mentum algebra [20].[21j. 

The complete set of eigenvalues of a SO a (3) multiplct 
is given by (|4"Tj) and (|42|) . Hence we are able to search 
for a realisation of SO a (3) symmetry in nature. Promis- 
ing candidates for a still unrevealed SO a (3) symmetry are 
the ground state band spectra of even-even nuclei. There- 
fore we will introduce the fractional symmetric rigid rotor 
model in the next section and discuss some of its applica- 
tions. 
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4 Eigenvalues of the fractional symmetric 
rigid rotor model 

The fractional Schrodinger equation is given by (see [7]) 

H& = E& (45) 

The requirement of invariance under fractional rotations 
in B? completely determines the structure of the fractional 
Hamiltonian up to a constant factor Aq: 

77 = m + A Q A\ (SO a (3)) (46) 
= m + A J 2 (a) (47) 

where mo mainly acts as a counter term for the zero point 
energy contribution of the fractional rotational energy and 
Aq is a measure for the level spacing. 

On a function set | JM > the Hamiltonian is diagonal. 
Furthermore, due to definitions the Hamiltonian 

commutes with the parity operator 77, [if, II] = 0. 

For a — 1 the function set | JM > reduces to the set 
of spherical harmonics Yjm in carthesian representation. 
The eigenvalues of the parity operator 77 are given by [27] 

nYjM(xi,X2,x 3 ) = Yjm(-x 1 ,-x 2 ,-x 3 ) 

= (~l) J Y JM (x 1 ,x 2 ,x 3 ) (48) 

Therefore, the wave function is invariant under parity trans- 
formation 77, if J is restricted to even, non-negative in- 
tegers J = 0, 2, 4, .... In a collective geometric model, this 
symmetry is interpreted as the geometry of the symmetric 
rigid rotor model [23] , 

Whether or not the behaviour ([48]) is still valid for the 
function set | JM > with arbitrary a cannot be proven 
directly with the methods developed so far. 

Nevertheless, restricting J to be an even, non-negative 
integer J = 0, 2, 4, ... for arbitrary a, implies a symmetry, 
which we call in analogy to the case a = 1 the fractional 
symmetric rigid rotor, even though this term lacks a direct 
geometric interpretation for a/1. 

Hence we define the spectrum of a fractional symmetric 
rigid rotor: 

r(l + {J+l)a) 

E = E°j =m +A h 2 — ) £ , 7 = 0,2,4,... 

r(l + (J-l)a) 

(49) 

In figure [1] this energy spectrum is plotted for different 
values of a. 

As a general trend the higher angular momentum en- 
ergy values are decreasing for a < 1. This behaviour is a 
well established observation for nuclear low energy rota- 
tional band structures [2"5] , 

In a classical geometric picture of the nucleus [22] , [23] 
this phenomenon is interpreted as a change of the nuclear 
shape under rotations, causing an increasing moment of 
inertia. Now, with a fixed moment of inertia ( Aq is a 
constant) in our approach for the fractional rotational en- 
ergy, the same effect results as an inherent property of the 
fractional derivatives angular momentum. 




Fig. 1. energy spectrum E J from (|49[) of the fractional sym- 
metric rigid rotor 



Another remarkable characteristic of the fractional sym- 
metric rigid rotor spectrum is due to the fact, that the 
relative spacing between different levels is changing (see 
figured]). 

In classical geometric models there are distinct analyt- 
ically solvable limits, e.g. the rotational, vibrational and 
so called 7-unstable limits [23] . 

These limits are characterised by their independence 
of the potential energy from 7. In that case, the five di- 
mensional Bohr Hamiltonian is separable [22], [26]. With 
the collective coordinates (3, 7 and the three Euler angles 
0i the product ansatz for the wave function 

*{p,iA) = fiP)*{iA) (so) 

leads to a differential equation for /3, expressed in canon- 
ical form 

E(n p ,T) (/3 2 /(/?)) 

The ground state spectra of nuclei, based on 7 indepen- 
dent potentials, are then determined by the conditions 

n p = Q (52) 
2r = J , r = 0,1,2, ... 

In the following we will prove, that the above mentioned 
limits are included within the fractional symmetric rigid 
rotor spectrum as special cases at distinct a- values. This 
is a further indication, that the fractional symmetric rigid 
rotor model may be successfully applied to low energy 
excitation spectra of nuclei. 
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Fig. 2. relative energy levels g(J, a) = {E a , - E%)/{E% - Eft) 
with EJ from (I49|l of the fractional symmetric rigid rotor. The 
thick line shows the minimum achievable ratio for a given J as 
a function of a. 



4.1 Rotational limit 

In the geometric collective model the rotational limit is 
described by the symmetric rigid rotor [24].|23|: 



J = 0,2,4,. 



(53) 



This limit is trivially included in the fractional symmetric 
rigid rotor spectrum for a = 1 



EJ =1 =m Q + A h 2 



r(i + (,/ + !)) 



mo + A Ti 



r(i + (j-i) 

2 (J + 1)! 
(J-l)! 



(54) 



= m + A h 2 J(J+ 1) , J = 0,2, 4, ... 

Setting mo = and A = 1/(20) completes the deriva- 
tion. 

Remarkably enough, reading (|54[) backwards is an ap- 
plication of the original Euler concept for a fractional 
derivative, since Euler suggested a replacement of n! for 
the discrete case by r(l + z) for the fractional case |17j . 



4.2 Vibrational limit 

In the geometric collective model the vibrational limit is 
described by a harmonic oscillator potential e. g. [23] 



(55) 



where C is the stiffness in /3-direction. The level spectrum 
is given by 

E = Uu}(N + 5/2) (56) 



with uj = y/C/B and 
N = 2np + t 



N = 0,1,2,3,. 



(57) 



Therefore the ground state band E 9 S - is given according 
to the conditions |52|) by 



E 9 - s -(t) = hu(r + 5/2) 



t = 0,1,2,. 



(58) 



We will prove, that for a = 1/2 the spectrum for the 
fractional symmetric rigid rotor corresponds to this vi- 
brational type ground state spectrum 



E 



a=l/2 
J 



mo + A h 



Wi + ±(j + i 



r i + i(J-r 



(59) 



= too + A h 



to + A h' 



2 r(i + j/2 + 1) 
r(i + j/2-i) 
2 r(i + j/2 + 1) 
r(j/2 + i) 



m + A r(J/2+ -) 



J = 0,2, 4, 



where we have used -T(l + z) = zT(z). 

Since J = 2r for the ground state band we get: 



a=l/2 
J 



1. 



to +A Q h 2 {T + -) 



0,1,2,... (60) 



Setting to = 2hu> and A = w/Ti completes the derivation. 

The appearance of the correct equidistant level spac- 
ing including the zero point energy contribution of the 
harmonic oscillator eigenvalues is due to the fact, that the 
Riemann fractional derivative does not vanish when ap- 
plied to a constant function. 

Instead, for the fractional symmetric rigid rotor based 
on the Riemann fractional derivative there always exists 
a zero point energy for a / 1 of the form 



i 2 Hloo>=^^ii-4 

r(l-a) 



00 > 



(61) 



Hence the consistence with the spectrum of the harmonic 
oscillator including the zero point energy is a strong ar- 
gument for our specific choice of the Riemann fractional 
derivative definition. 



4.3 Davidson potential - the so called -unstable limit 

There exists a third analytically solvable case, which is 
based on the Davidson potential [25] , which originally was 
proposed to describe an interaction in diatomic molecules. 
The potential is of the form 



™-s c *(|-f 



(62) 



where (3q (the position of the minimum) and C (the stiff- 
ness in /3-direction at the minimum) are the parameters 
of the model. 
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For Pq = this potential is equivalent to the harmonic 
oscillator potential. 

The energy level spectrum is given by 



E{n ,r) = hu{n + - + -a T ) - -Cffi 



with 



13CI% + (2t + 3) 2 



For the ground state band we get 



1 1 



E 3 ' s '(t) = fuj[- + -a T — 7CA 



(63) 
(64) 

(65) 



For > we expand the square root in a T in a Taylor- 
Series 



Shifting t by 



= f - 3/2 



(66) 



(67) 



causes the linear term to vanish and the resulting level 
scheme is of the form 



4.4 Linear potential limit 

Besides the above discussed well established three limits 
there exists an additional special case, which is, in a geo- 
metric picture, the linear potential model. It has first been 
proposed by Fortunato|26j. 

Although it has not been used for a description of nu- 
clear ground state band spectra yet, it will turn out to be 
useful for an understanding of nuclear spectra near closed 
shells discussed in the following section. 

The potential is of the form 



V(p) = C/3 



(71) 



where C (the stiffness in /3-direction) is the main param- 
eter of the model. 

The level spectrum is given approximately by 26 : 



3 1 /3C 

E{np,T) = -C/3 + (n /3 + -)W — 



with 



0o= (2(t + 1)(t + 2)/<7) 



1/3 



(72) 



(73) 



The ground state band level spectrum is given according 
to conditions (|52|) by 



E 9 - s -(t) 



Cq + CiT +.. 



(68) 



Therefore the 7-unstable Davidson potential is charac- 
terised by the condition, that the linear term in a shifted 
series expansion in r vanishes. 

In order to determine the corresponding fractional co- 
efficient a, we shift the fractional energy spectrum by 
—3/2 and expand in a Taylor series at J/2 = r = 
and denote the di- and trigamma function): 



, E(l - a/2) 
r 2 1 (1 — 5a/2) 

1 - a ($(1 - ha/2) - V(l - a/2) J r 
^{((V(l-5a/2)-V(l-a/2)) 



(69) 



W 1 ^ - ha/2) -W 1 ^- a/2) 



>T* + .. 



The linear term in r has to vanish. Therefore a is deter- 
mined by the condition 



&(l-ha/2) =#(1 -a/2) 



(70) 



which is fulfilled for a ~ 0.66. 

Hence for a sufficiently large /3o the ground state band 
spectrum of the 7-unstable Davidson potential is repro- 
duced within the fractional symmetric rigid rotor at a ~ 
2/3. 



(74) 



We therefore are able to define the relative energy levels 
/(r) in units E g s (l) — E 9 S (0), which are independent of 
parameter C. 



fir) 



E3- s -(t) - E9- s -(0) 



(75) 



E9- s -(l) - E9- s -(0) 

An expansion in a Taylor series at r — 1 yields: 

f(r) = 1 + 0.895(t-1)-0.076(t-1) 2 + 0.018(t-1) 3 + ... 

(76) 

An equivalent expression for the fractional rotational en- 
ergy is given by: 



50, a) 



E% 



Eg 



(77) 



A Taylor series expansion at r = 1 followed by a compar- 
ison of the linear terms of /(r) and g(r, a) leads to the 
condition 



2<xT(l - a)r(3a) (2 + 3a<P{a) - &(3a) 
r 2 {l + a) - T(l - a)r(l + 3a) 



0.895 (78) 



which is fulfilled for a 0.33. Hence below the vibrational 
region at a « 1/2 there exists a region at a w 1/3, which 
corresponds to the linear potential model in a geometric 
picture. 

Summarising the results of this section, we conclude, 
that a change of the fractional derivative coefficient a may 
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be interpreted within a geometric collective model as a 
change of the potential energy surface. 

In a generalised, unique approach the fractional sym- 
metric rigid rotor treats rotations at a » 1, the 7-unstable 
limit at a ~ 2/3, vibrations at a « 1/2, and the linear po- 
tential limit at a w 1/3 similarly as fractional rotations. 
They all are included in the same symmetry group, the 
fractional SO a (3). This is an encouraging unifying point 
of view and a new powerful approach for the interpretation 
of nuclear ground state band spectra. 



5 Comparison with experimental data 

In the previous section we have shown, that the rotational, 
vibrational and 7-unstable limit of geometric collective 
models are special cases of the fractional symmetric rigid 
rotor spectrum. 

The fractional derivative coefficient a acts like an order 
parameter and allows a smooth transition between these 
idealised cases. 

Of course, a smooth transition between rotational and 
e.g. vibrational spectra may be achieved by geometric 
collective models too. A typical example is the Gneuss- 
Greiner model [28] with a more sophisticated potential. 
Critical phase transitions from vibrational to rotational 
states have been studied for decades using e.g. coherent 
states formalism or within the framework of the IBA- 
model [29]- [32]. But in general, within these models, re- 
sults may be obtained only numerically with extensive ef- 
fort. 

Anyhow, only very few nuclear spectra can be described 
accurately by one of the four limiting, idealised cases dis- 
cussed in section [4] 

In the following we will prove, that the full range of 
low energy ground state band spectra of even-even nu- 
clei is accurately reproduced within the framework of the 
fractional symmetric rigid rotor model. 

For a fit of the experimental spectra Ej Xp we will use 



Ej — mo + a 



r(l + (J+l)a 
r(i + (J - l)a 



J = 0,2,4,.. 



(79) 

with the slight modification, that fi has been included in 
to the definition of ao, so that m and a may both be 
given in units of [/ceV]. 

As a first application we will analyse typical rotational, 
7-unstable, vibrational and linear type spectra. In the up- 
per row of figure [3] the energy levels of the ground state 
bands are plotted for 156 Gd, 196 Pt, 110 Cd and 218 Po which 
represent typical rotational-, 7-unstable-, vibrational and 
linear type spectra [33]-|36j. 

The fractional coefficients a, deduced from the experi- 
mental data, are remarkably close to the theoretically ex- 
pected idealised limits of the fractional symmetric rigid 
rotor, discussed in the previous section. 



Table 1. Listed are the optimum parameter sets (a, ao, mo 
according (|79|> ) for the fractional symmetric rigid rotor for dif- 
ferent nuclids. The maximum valid angular momentum J max 
below the onset of alignment effects are given as well as the 
root mean square error AE between experimental f |33J-[55] ) 
and fitted energies in %. 



nuclid 



a [keV] m [keV] J max AE[%] 



92 Gde4 


0.863 


31.90 


-43.65 


14 


2.23 


196 p. 
118^*78 


0.710 


175.06 


-83.69 


10 


0.44 


1 62Cd48 


0.570 


607.91 


-405.80 


6 





218p„ 
134^0-84 


0.345 


1035.69 


-671.03 


8 


0.12 


"""llOsre 


0.624 


339.423 


-128.448 


6 





^Osre 


0.743 


125.128 


-39.968 


10 


1.35 


960S76 


0.767 


89.010 


-23.656 


8 


1.21 


980S76 


0.771 


63.388 


-17.656 


24 


0.73 


10oOs76 


0.808 


51.231 


-23.064 


18 


1.22 


lolOs76 


0.816 


49.882 


-22.792 


14 


2.16 


180 r>o 
IO4US76 


0.841 


45.309 


-18.662 


12 


2.50 


106OS76 


0.903 


32.002 


-7.159 


10 


1.29 


108OS76 


0.904 


31.690 


-12.902 


14 


1.39 


II0OS76 


0.882 


40.433 


-19.155 


14 


1.82 


II2OS76 


0.875 


44.567 


-14.629 


12 


1.58 


190 n c , 
114 US7 6 


0.847 


58.055 


-17.748 


12 


1.34 


116 ( JS70 


0.835 


63.774 


-15.437 


12 


0.71 



_>14 

120 



120 



Rass 



-0.007 
0.548 



374408 
344.107 



-376529 
305.766 



8 
24 



2.53 
8.27 



?2sRa88 

?3oR- a 88 
220 t 
132-1 
222 1 
134-1 

13 6 Ra88 



0.181 
0.536 
0.696 
0.831 
0.841 



3665.36 
321.622 
83.603 
33.221 
27.295 



-2887.88 
-160.787 
-25.966 
-5.67 
-8.849 



10 
30 
30 
6 

12 



4.22 
1.25 
0.26 


1.57 



Table 2. Energy levels for 176 0s with optimum parameter set 
from table[T]for the fractional symmetric rigid rotor EJ accord- 



ing to (|79[) compared with the experimental values Ej from 
[41] and relative error AE in % for different angular momenta 
J. Note that for J > 147! beginning microscopic alignment ef- 
fects are ignored in the fractional symmetric rigid rotor model 
and therefore the error is increasing. 



j[h] 


Ej[keV] 


Ef p [keV\ 


AE[%] 


0+ 


-13.1 


0.0 




2+ 


144.9 


135.1 


-7.29 


4+ 


404.4 


395.3 


-2.31 


6+ 


744.8 


742.3 


0.33 


8+ 


1155.4 


1157.5 


0.17 


10+ 


1629.6 


1633.8 


0.25 


12+ 


2162.3 


2167.7 


0.24 


14+ 


2749.8 


2754.6 


0.17 


(16+) 


3389.2 


3381.4 


-0.23 


(18+) 


4077.8 


4019.1 


-1.46 


(20+) 


4813.6 


4683.2 


-2.78 


(22+) 


5594.8 


5398.8 


-3.63 
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Fig. 3. The upper row shows energy levels of the ground state bands for Gd, Pt, Cd and Po are plotted for 
increasing J. Squares indicate the experimental values given in [33 - 36 , +-symbols indicate the optimum fit of the classical 
symmetric rigid rotor and the curves give the best fit for the fractional symmetric rigid rotor. In the lower row the corresponding 
backbending plots are shown. 



From the experimental data, we can roughly distin- 
guish a rotational region for 1 > a > 0.8, a 7-unstable 
region for 0.8 > a > 0.6, a vibrational region for 0.6 > 
a > 0.4 and a linear potential region for 0.4 > a > 0.2. 
Experimental results are reproduced very well for all re- 
gions within the framework of the fractional symmetric 
rigid rotor model, while a similar fit with the classical 
symmetric rigid rotor gives much worse results, especially 
for the vibrator-like and linear spectra. 

The lowering of higher angular momentum states is 
correctly described for rotational-, 7-unstable-, vibrational- 
as well as linear type spectra by the fractional symmetric 
rigid rotor model. 

Since for higher angular momenta commencing micro- 
scopic effects limit the validity of the macroscopic frac- 
tional symmetric rigid rotor model fits, in the lower row 
of figure[3]the corresponding backbending plots are shown. 
With these plots, the maximum angular momentum J max 
for a valid fit may be determined. In table [1] the optimum 
parameter sets (ev, ao, mo) as well as the average root mean 
square deviation AE are listed. In general, the difference 
between experimental and calculated energies is less than 
2%. 

This means, that the error is better than a factor 3 
to 6 compared to the results of a standard Taylor series 
expansion up to second order in J. 

As a second application of the fractional symmetric 
rigid rotor, we study systematic isotopic effects for osmium 
isotopes (37]- [49]. Optimum fit parameter sets are given in 
the second part of table [H the corresponding backbending 
plots are given in figure [4] In table [2] the results for the 
fitted ground state band energies of the fractional sym- 
metric rigid rotor model according (|79[) are compared to 
the experimental levels [41] for 176 0s. 

Obviously below the beginning of alignment effects 
the spectra are reproduced very well within the fractional 
symmetric rigid rotor model. The fractional coefficient a 
listed in table [1] increases slowly for increasing nuclcon 
numbers. Consequently, we observe a smooth transition 



of the osmium isotopes from the 7-unstable to the rota- 
tional region. 

As a third application, we investigate the systemat- 
ics of the ground state energy level structure near closed 
shells. For that purpose, in figure [5] a series of spectra of 
radium isotopes near the magic neutron number N = 126 
are plotted, experimental data are taken from |50]-[55|. 

Starting with 224 Ra, the full variety of possible spectral 
types emerges while approaching the magic 214 Ra. 

While 224 Ra and 222 Ra are pure rotors, 220 Ra shows a 
perfect 7-unstable type spectrum. 218 Ra presents a spec- 
trum of a almost ideal vibrational type. 216 Ra is even 
closer to the magic 214 Ra and represents a new class of 
nuclear spectra. 

In a geometric picture, this spectrum may be best in- 
terpreted as a linear potential spectrum as proposed in 
section l4~4l Typical candidates for this kind of spectrum 
are nuclei in close vicinity of magic numbers like 2 ||Pog4, 

154VK 134v„ 967 v 88c,, 
84 YD 70, 80 Ae 54, 56^40 Or 50 Ol 38 . 

Finally, the experimental ground state spectrum of the 
magic nucleus 214 Ra shows a clustering of energy values: 
For low angular momenta a tends towards zero, which 
becomes manifest through almost degenerated energy lev- 
els while for higher angular momenta the spectrum tends 
to the vibrational type. In a classical picture this is in- 
terpreted as the dominating influence of microscopic shell 
effects. Within the framework of catastrophe theory [SB] 
this observation could be interpreted as a bifurcation as 
well. 

Therefore the fractional symmetric rigid rotor model 
is not well suited for a description of the full spectrum of 
a nucleus with magic proton or neutron numbers. But all 
other nuclear spectra are described with a high grade of 
accuracy within the framework of the fractional symmetric 
rigid rotor. 

As a remarkable fact a reduces smoothly within the 
interval 1 > a > while approaching a magic number, e.g. 
N = 126. Therefore, a is an appropriate order parameter 
for such sequences of ground state band spectra. 
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0.00 C.D2 MM 0.06 08 0.10 0.12 O il 0.00 0.05 0.10 0.15 20 0.25 0.30 0.00 0.D5 0.10 0.15 20 0.25 0.30 0.35 




0.01 0.02 0.03 0.01 0.05 0.06 0.07 0.08 0.09 0.10 0.11 qui 002 0.03 O.Oi D.05 0.0B 0.07 O.QS 0.09 0.10 OM 0.12 0.13 0-0' D.02 0.03 D.Oi 0.05 0.0B 0.07 0.0B 0.09 0.10 0.11 0.1! 0.13 



Fig. 4. For a set of osmium isotopes backbending plots are shown. Squares indicate experimental values, given in [37|-|49], 
The thick line is the optimum fit result for the fractional symmetric rigid rotor. Optimum fit parameter sets are given in table 
[1] For the fit all angular momenta below the onset of microscopic alignment effects are included. 




10 20 j 30 2 i I I j 10 2 * 6 8 10 12 U 16 18 20 22 2( 



Fig. 5. The change of energy level structure of ground state bands near closed shells is illustrated for the magic neutron number 
N = 126 with a set of radium isotopes. Squares indicate experimental energy values, given in [50]-[55]. Crosses indicate the 
optimum fit with the standard symmetric rigid rotor model. The thick line is the optimum fit result for the fractional symmetric 
rigid rotor. Optimum fit parameter sets are given in the lower part of table [T] For 214 Ra, two different fits for low and higher 
angular momentum indicate the different spectral regions for magic nucleon numbers. 
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100 



80 



60 



20 



0.0 0.1 2 0.3 0.4 0.5 0.6 7 0.8 0.9 1.0 

a 

Fig. 6. The distribution n(a) of a-values from a fit of 490 
ground state band spectra of even-even nuclei, accumulated 
in intervals Act = 0.1 from |Be4 to f^gFmioo and 152N0102 
respectively. 




Fig. 7. Plot of fitted a parameters for a set of osmium isotopes 
as a function of increasing neutron number N. Squares indicate 
the optimum a listed in table [1] Thick lines are fits according 
to equation J5T]). 



As a fourth application we analyse some general as- 
pects, which are common to all even-even nuclei. Experi- 
mental ground state band spectra up to at least J = 4+ 
are currently known for 490 even-even nuclei [5 7) . For 
these nuclei a full parameter fit according ([79]) was per- 
formed. The resulting distribution n(a) of fitted a values 
in 0.1 steps is plotted in figure [TJ] 

This distribution is not evenly spread. Most of even- 
even nuclei exhibit rotational type ground state band spec- 
tra followed by 7-unstable and vibrational type spectra. 
There is a distribution gap at a ~ 0.3. Remarkably enough, 
this is the region, where in a geometrical picture nuclear 
spectra corresponding to the linear potential model are 
expected. 

About 8% of even-even nuclear spectra are best de- 
scribed with a w 0.2. 

The sequence of radium isotopes, presented in figure 
[5] reflects the general distribution of spectra for even-even 
nuclei. 

An interesting feature may be deduced from the ob- 
servation, that n(a) is almost linear in the region 0.3 < 
a < 0.8. 

Hence, with the proton number Z and the neutron 
number N, we define the distance R(Z, N) in the (Z, N)- 
plane, from the next magic proton or neutron number re- 
spectively to be 



R(Z, N) = y '{Z - Z mag f + (N - N mag f 



(80) 



The linearity of n(a) implies a linear dependency of a on 
R{Z,N): 

a(Z,N) = c + c 1 R(Z,N) , 0.25 < a < 0.85 (81) 

which is a helpful relation to determine an estimate of a 
for a series of nuclei. 

Therefore the series of a values for osmium isotopes, 
given m tabled] may be understood even quantitatively: 



The series of Os to Jj^Os is closer to the N = 82 
magic shell and therefore shows an increasing sequence in 
a. A least square fit yields c — 0.528 and ci = 0.0143 
with an average error in a less than 3% 

The series JgfOs to }?§Os is closer to the N = 126 shell 
and therefore shows a descending sequence in a. A least 
square fit yields Co = 0.744 and c\ — 0.008 with an error 
in a less than 1%. 

In figure[7]these results are presented graphically. From 
the figure we can deduce a prediction for cv^gpOs) = 0.671 
and a( gfOs) = 0.695. From these predicted a values the 
ratio 



TPOt 

h 4+ 



rpry 



Tp a. rp a 



(82) 



may be deduced. 



For §gOs a value of R^ + /2+ 



2.29 and for ^fOs a 
value of i? 4+ / 2 + = 2.35 results. These predictions are yet 
to be verified by experiment . 

Furthermore, the sequence of osmium isotopes, pre- 
sented in figure [7] indicates the position of the next magic 
neutron number N — 126. The nucleus If&Os with the 
maximum a value in the sequence of osmium isotopes, is 
positioned just between the two magic numbers N = 82 
and N = 126 and may be used for an estimate for the 
next magic neutron number. 

Therefore, we propose an alternative approach to the 
search for super heavy elements. Instead of a direct syn- 
thesis of superheavy elements [58], [59], a systematic sur- 
vey of a values of the ground state band spectra of even- 
even nuclei with N > 146 and Z > 98 could reveal indi- 
rect information about the next expected magic numbers 
Z = 114 and N — 164. The currently available experimen- 
tal data do not suffice to determine these magic numbers 
accurately. 

From the results presented we draw the conclusion, 
that the full variety of low energy ground state spectra of 
even-even nuclei is described with a high grade of accuracy 
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from a generalised point of view within the framework of 
the fractional symmetric rigid rotor. 

The great advantage of the fractional symmetric rigid 
rotor model compared to the classical geometric models 
is, that nuclear ground state band spectra are described 
analytically with minimal effort according to equation (|T9[) 
and with an excellent accuracy. 

6 Conclusion 

The Casimir operators and multiplets for a fractional ex- 
tension of the rotation group SO(n) have been classified 
algebraically. 

The major result of a discussion of the corresponding 
fractional symmetric rigid rotor spectrum is, that the ro- 
tational, vibrational and so called 7-unstable limit of geo- 
metric collective models are special cases of the fractional 
symmetric rigid rotor spectrum. They are all treated as 
generalised rotations and are included in the same sym- 
metry group, the fractional SO a (3). 

The fractional derivative coefficient a acts like an order 
parameter and allows smooth transitions between these 
idealised cases. 

Applied to the ground state band spectra of nuclei, 
there is an excellent agreement with the experimental data 
for the full range of angular momenta up to the begin- 
ning of alignment effects. The results indicate, that the 
fractional symmetric rigid rotor is an appropriate tool to 
reproduce the low energy ground state excitation spectra 
of even-even nuclei. 

The results encourage further studies with this model. 
With the knowledge of the corresponding eigenfunctions 
of the fractional symmetric rigid rotor e.g. a calculation 
of transition matrix elements was possible. 
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